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ABSTRACT 

V893 Sco is an eclipsing dwarf nova that had attracted little attention from 
X-ray astronomers until it was proposed as the identification of an RXTE all-sky 
slew survey (XSS) source. Here we report on the pointed X-ray observations of 
this object using Suzaku. V893 Sco was in quiescence at the time, as indicated 
by the coordinated optical photometry we obtained at the South African Astro- 
nomical Observatory. Our Suzaku data show V893 Sco to be X-ray bright, with 
a highly absorbed spectrum. Most importantly, we have discovered a partial 
X-ray eclipse in V893 Sco. This is the first time that a partial eclipse is seen 
in X-ray light curves of a dwarf nova. Our preliminary simulations demonstrate 
that the partial X-ray eclipse can be in principle reproduced if the white dwarf 
in V893 Sco is partially eclipsed. Higher quality observations of this object have 
the potential to place significant constraints on the latitudinal extent of the X- 
ray emission region and thereby discriminating between an equatorial boundary 
layer and a spherical corona. The partial X-ray eclipse therefore makes V893 Sco 
a key object in understanding the physics of accretion in quiescent dwarf nova. 
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Subject headings: Stars: novae, cataclysmic variables — stars: individual (V893 
Sco) — X-rays: binaries 



Introduction 



Cataclysmic variables (CVs ), in which a w hite dwarf primary accretes from a Roche-lobe 
filling, late- type secondary (see IWarner Ill995l for a comprehensive review), are an excellent 
laboratory for the physics of accretion. In the subclass of dwarf novae, accretion proceeds 
via a disk that switches between the low (quiescence) and high (outburst) luminosity states. 
The visible light of a quiescent dwarf nova is usually dominated by the bright spot, where the 
accretion flow from the secondary hits the outer edge of the disk, as well as the photosphere 
of the white dwarf. In outburst, the disk becomes the dominant source of visible light. 
In contrast. X-ray observations of dwarf novae probe the accretio n flow in the immediate 



vicinity of the white dwarf, such as an optically thin boundary layer (jPatterson fc Raymond 



19851 ). since the accretion disk in a dwarf nova is too cool to emit X-rays. Although the disk 
instability model is a highly successful framework for explaining the dwarf r iova outbursts , 
there are details that defy the prediction of the basic version of the model (ILasota II2OOII ). 
In partic ular, the observed X -ray luminosities of quiescent dwarf nova (often of order 10^^ 



ergss"^; iBaskill et al. I l2005l ) imply an accretion rate onto the white dwarf that is much 
higher than predict ed. Proposed modifications of th e disk instability model include the 
coronal siphon flow (IMeyer fc Meyer- Hofmeister Ill994j ). which might lead to accretion over 
a much wider area of t he white dwarf surface than through a boundary layer, and a weakly 
magnetic white dwarf (ILivio fc Pringle Ill992i ). In the latter model, the magnetic field is 
too weak to control the accretion flow during outburst, but strong enough to do so during 
quiescence. 



Kato et al. 



V893 Sco is an eclipsing (IBruch et al. Il2000l) dwa rf nova with an orbital period of 1.8 hr 



( iThorstensen Ill999l ). According to 
of ~14.5, and has an outburst every 
Exte nsive photometry by Warner and collaborators (j Warner et al. I l2003l : iPretorius et al. 



2002), V893 Sco has a quiescent magnitude 
30 days duri ng which it reaches magnitude ~12. 5. 



20061 ) have revealed quasi-periodic oscillations and dwarf nova oscillations in V893 Sco, but 
a strictly periodic signal was never found. Therefore, existing optical data point strongly 
towards a dwa rf nova classifica tion and argue against an intermediate polar (IP, or DQ Her 
type systems; [Patterson 1 119941 ) classification. IPs are a subset of magnetic CVs in which 
the primary's magnetic field disrupts the inner accretion disk, channeling the flow to the 
magnetic polar region(s); the spin period of the magnetic white dwarf is a strict clock that 
characterize IPs. Thus, V893 Sco joins a growing number of eclipsing dwarf novae below the 
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period gap. 



Citing the vari able shape of the ecli pses and the relatively small eclipse amphtude (often 
less than 0.75 mag). lBruch et al. I (120001 ) argue for a gr azing eclipse of the bright s pot and the 
disk, but not of the white dwarf. On the other hand, iMatsumoto et al. I (120001 ) argue that 
the white dwarf is eclipsed, based on the fact that the spectroscopic conjunction of the disk 
(the red-to-blue crossing of the emission line radial velocities) occurs at mid-eclipse. From 
the radial velocity curves, they also infer a wh ite dwarf mass of ~0.5-0.6 Mq and a mass 
ratio of ~0.2-0.3. However, iMason et al. I (120011 ) prefer a higher mass (0.89 M©) white dwarf 
and an inclination angle of 72.5 degrees based on their analysis of the emission line radial 
velocities, although this value depends in part on the assumed mass-radius relationship for 
the secondary. 

Despite the current lack of a consensus regarding the system parameters, the eclipsing 
nature of V893 Sco makes it an important target for detail ed studies. Moreov er, it is a 



nearby system with parallax-estimated distance of 155_34 pc (IThorstensen II2003I ). 



Although the ROSAT detection was already noted by lKato et al. I (119981 ) . V893 Sco did 
not draw the attentio n of X-ray astronorners un til it was listed among the RXTE all-sky slew 
survey (XSS) sources (IRevnivtsev et al. II2004I ). along with three other dwarf novae (SS Aur, 
V426 Oph, and SU UMa). The estimated luminosities of these 4 systems are all just under 
Xq32 gj-ggg-i in the 2-10 keV band, placing them near the upper end for no n-magnetic CVs 



10^°-3 X 10^2 in 0.1-100 keV for the ASCA sample; iBaskill et al. 1120051 ). However, the 



XSS is based on the data taken during slews of the non-imaging RXTE PCA detector. For 
relatively faint sources such as these dwarf novae, the positional errors are considerable, and 
misidentification is a possibility. Hence it is important to confirm the proposed identification 
of XSS dwarf novae. This is particularly true for V893 Sco, which had never been the subject 
of a pointed X-ray observation above 2 keV. Although V893 Sco is securely detected as a 
ROSAT all-sky survey (RASS) source, the ratio of XSS (2.38 c/s in the 3-8 keV band) 
to RASS (0.35) count rates is high, compared to, e.g., SS Aur which has the same RASS 
count rate but was detected at 0.75 c/s in the XSS 3-8 keV band. This leaves open the 
possibility that V893 Sco is only a partial contributor to the XSS flux. For this reason and 
also because of the potential return in studying the X-ray properties of this eclipsing dwarf 
nova, we performed a pointed observation using Suzaku, along with contemporaneous optical 
photometry. 
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Observations and Data Reduction 



We observed V893 Sco with Suzaku flMitsuda et al. Il2007f ) between 2006 Aug 26 09:40 



and Aug 27 01:20 UT (sequence number 401041010). There are two types of co-aligned 
instruments on-board Suzaku, which are: 4 units of the CCD-based X-ray imaging spec- 
trometer (XIS), each behind its X-ray telescope (XRT), and the non-imaging hard X-ray 
detector (HXD). As the source is relatively faint for a non-imaging instrument and the ex- 
posure time short, we did not analyze the HXD data, a decision that we justify at the end 
of §3.1. 

We started from data processed using the V2.0.6.13 pipeline. We then updated the 
energy scale calibration using the 2008 Feb 1 release of the calibration database. Initially we 
applied the standard screening criteria: attitude control system in nominal mode, pointing 
within 1.5 arcmin of the mean direction, XIS data rate medium or high, the satellite outside 
the South Atlantic Anomaly (SAA), at least 436 s after the last SAA passage, elevation above 
Earth limb > 5°, elevation about the bright Earth limb > 20°. We then experimented with 
relaxing some of the standard criteria to maximize the phase coverage of the observation. We 
found that relaxing the post-SAA passage limit from 436 s to 180 s, and the bright Earth 
limb limit from 20° to 15° did not noticeably increase the noise level of our data. With 
this modification, the total XIS exposure time was 19,174 s, representing a 670 s increase 
compared to the standard screening. The background level during the screened interval was 
normal for Suzaku, and varied within ±50% of the mean value primarily as a function of the 
geomagnetic coordinates of the spacecraft. 

We found a bright source centered 13.5 arcsec from the precisely known optical position 
of V893 Sco in the XIS image. Since the current r econstruction of Suzaku attitude results 



in a 90% error circle of 19 arcsec for bright sources (jUchiyama et al. II2008I ). we identify the 
Suzaku source with V893 Sco. Given the statistical quality of the image, the apparent offset 
is caused predominantly by the attitude reconstruction errors. 

We used a 3.5 arcmin radius extraction region for the source and 7.5 and 4 arcmin radius 
annular extraction region for the background, both centered on the location of V893 Sco as 
measured on the XIS image. For spectroscopy, we summed the data and the responses of 
the three XIS units with frontside illuminated (FI) CCD chips, because they have nearly 
identical responses. The mean net rate of V893 Sco was 0.85 cts~^ per FI unit, while the 
estimated background rate was 0.04 cts~^ in the source region. The data from XISl (with 
the backside illuminated, or BI, CCD chip) have a higher background rate averaging 0.16 
cts~^, although the net source rate was also higher at 0.98 cts~^. We fit the XISl spectrum 
simultaneously with the sum of FI spectra, after grouping channels so they have a minimum 
of 25 counts, thus making the statistic an acceptable approximation. For photometry. 



- 5 - 



we added background subtracted light curves from all 4 XIS units over the 0.4 keV (FI)/0.3 
keV (BI) to 10 keV bandpass using 8 s bins, preserving the original time resolution of the 
XIS data. We also extracted light curves below 2 keV, in the 2-4 keV band, and above 4 
keV to investigate possible energy dependence. 

We also carried out optical photometry of V893 Sco on the nights of 2 5, 26, and 27 



Augu st 2006 using the University of Cape Town CCD Photometer (UCTCCD; lO'Donoghue 



19951 ) on the South African Astronomical Observatory (SAAO) 1.9m telescope. The detec- 
tor is a Wright Instruments Peltier-cooled CCD with a thinned and back-illuminated EEV 
P86321/T chip. It was used in frame-transfer mode where one half of the CCD is masked 
and used as a readout area, allowing exposures to be made with essentially no "dead" time 
between exposures. On the 1.9m telescope, the pixels on the detector correspond to 0.13 
arcsec so that it is normal to use 3x3 or greater pre-binning, unless the seeing is better than 
about 1 arcsec, to ensure optimal data-extraction. Because V893 Sco is relatively bright, 
integration times of 6 seconds were used on all nights. All observations were made in white 
light. The CCD was pre-binned at 5x5, 4x4, and 3x3 for the three nights respectively. 

Observations were reduced at the telescope using a custom pipehne, which allows the 
observer to select suitable comparison stars and judge observations. The conventional proce- 
dure of reducing observations were followed (bias subtraction, flatfield corr ection etc.) after 



which the stellar brightnesses were extracted using the method described in lSchechter et al 



(119931 ). Differential corrections were made by using the lightcurve of the brightest nearby 
star. Due to the small size of the UCTCCD (50x34 arcsec^ on the 1.9m telescope) the only 
available comparison stars were somewhat fainter than V893 Sco but they were nevertheless 
used as the differentially corrected lightcurves were deemed to be of better quality than the 
uncorrected ones. With this set-up, the statistical fluctuation of the data is much less than 
the intrinsic variability of V893 Sco. Typical number of counts in a 6 s integration was 
100,000, and it never dropped below 30,000 even during the echpse. 



3. Results 



3.1. X-ray Spectrum 



Low resolution X-ray spectra of dwarf novae are generally dominated b y an absorbed 
brenisstrahlung-like conti nuum and a prominent He- like Fe Ka at 6.7 keV (IBaskill et al. 
20051 : iPandel et al. II2005I ). The Suzaku spectrum of V893 Sco is no exception, as shown in 
Figuredl The Fe Ka line indicates the origin in an optically-thin, thermal plasma. Therefore 
the simplest model consistent with the data is the single temperature plasma model mekal 
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( Mewe et al. 1985 . 1986 : Liedahl et al. Ill995 : Kaastra et al. 1996 ) as implemented in Xspec 
flArnaud Ill996f )7 There is a strong reason to believe that this is not a physically correct model 
for X-ray emission in CVs: the post-shock plasma in CVs loses energy by emitting X-rays 
and therefore must cool, and cannot remain at a single temperature. However, it is certainly 
possible for a multi-temperature plasma to have a strong peak in the differential emission 
measure distribution such that a single-temperature plasma model provides a good fit. For 
this reason, and to allow easy comparison with other works, we present the results of the 
mekal model fit in Tabled! 

However, as our primary spectral model, we choose t he mkcf low model that was orig i- 
nally developed for the cooling flow in clusters of galaxies (IMushotzkv fc Svmkowiak 1119881) . 
but was successfully applied to X-ray emission from CVs (IMukai et al. II2003I ). In its original 
context, the cooling flow model assumes a steady-state with a large reservoir of plasma in 
the outer parts of the cluster at the maximum temperature, kT^a^; the plasma cools by 
emitting X-rays with no additional energy input. This model is applicable to CVs as long 
as the post-shock region is in a steady state and there are no additional heating or cooling 
terms beyond X-ray emission. 

There are oth er models of multi-temperature X-ray emissi on specially developed for 



magne tic CVs (e.g.. lCropper et al. I I1999I : ISuleimanov et al. 1 l2005h . However, as lMushotzky &: Symkowiak 



( 119881 ) state (section HI): "Assuming that the same mass flow rate pertains throughout the 
cooling flow, the emission measure for each temperature is determined by the time it takes 
for the matter to radiate away sufficient energy to cool down to the next temperature shell. 
The differential emission measure is thus proportional to the reciprocal of the bolometric 
luminosity at that temperature." It is therefore unnecessary to solve for the structure of 
the emission region in order to predict the correct output X-ray spectrum, in the absence of 
additional heating or cooling terms. Thus, mkcf low serves well a 1st order approximation, 
both in clusters of galaxies and in CVs. 

We present the results of our fit in Figured] (for the mkcf low model) and in Tabled] (for 
both models). Note that we have used a global abundance value in our flts, and in effect it 
is the Fe abundance we have measured. This is because the only discrete feature evident in 
the Suzaku spectrum of V893 Sco is the Fe Ka complex; the abundances of other elements 
cannot be constrained reliably, in the absence of discrete features. To investigate the Fe 
Ka complex further, we fltted an absorbed bremsstrahlung continuum and three Gaussians 
in the 5-10 keV range. We do not detect any energy shifts or broadening of the Fe lines, 
although the upper limits are unconstraining due mainly to the modest signal-to-noise ratio 
of the XIS spectra. In addition to the dominant 6.7 keV (He-like) line with an equivalent 
width of 400±30 eV, we also detected both the H-like fine (6.97 keV; 100±20 eV) and the 
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fluorescent line at 6.4 keV (40±15 eV). The small equivalent width of the 6.4 keV line 
indicates that optically thick matter (white dwarf surface and an optically thick inner disk; 
Done &: Osborne 1 119971 ) covers less than Ivr steradian of the sky as seen from the emission 



region, which would result in an equivalent width in the 100-150 eV range (iGeorge fc Fabian 



19911 ). Although we did not include the fluorescent line in the global fit, it is sufficiently 
weak not to have affected the quality of the fit or the derived parameter values. 

We note that the derived abundances depend on our choice of spectral model (Tabled]). 
It is likely that a single temperature mekal model underestimates the Fe abundance, because 
the best fit temperature is near the peak of He-like Fe Ka emissivity. In the mkcf low fit, the 
plasma temperature range extends well outside the peak emissivity of He-like Fe Ka, hence 
a higher abundance is required to fit the observed data, and this probably is a more realistic 
refiection of the true abundances of V893 Sco. 

With both mekal and mkcf low models, the use of a single absorber (using the wabs 
model) resulted in significant residuals particularly below 1 keV (see the middle panel of 
Figured]). The fitted Nh of > 5 x 10^^ cm~^ is too high to be the interstellar column towards 
y893 Sco at d=155 pc; eve n the total Galactic column in this direction is only 1.4x10^^ cm~^ 
( Dickey &: Lockman 1990) as reported by the HEASARC Nh too|l|. Thus, we conclude that 
the measured column is dominated by an intrinsic absorber within the V893 Sco binary 
system. For e xample, our line of sight to the boundary layer may pass through the inner 



accretion disk (Ivan Teeseling et al. Ill996l ) 



The magnitude of the low energy residuals is far too large to be removed by changing the 
assumed abundances of the absorber. For example, substituting wabs wit h phabs model in- 
stead , and tr ying various "solar" abundance values (e.g., those according to I Anders fc Grevesse 
19891 and by IWilms et al. I l2000l ) with the mkcf low model resulted in essentially the same 
values of reduced (1.32-1.34). One can get an improved fit (reduced x^=l-09) by chang- 
ing abundances element by element using the varabs model, but at the expense of non- 
astrophysical abundances (>20 times solar for Na, for Al through Cr, and essentially solar 
for Fe). On the contrary, the low energy residual was successfully removed by adding an 
additional partial covering absorber (pcf abs) (see the top and the bottom panels of Figured] 
and Tabled])- 

We also tried using an additional emission component instead, while keeping just a 
single wabs absorber. These were not as successful at removing the residuals as the complex 
absorber: adding a low temperature mekal to mkcf low resulted in x^=1.08 with residuals 
below 1 keV, and the addition of a blackbody resulted in x^=1.04 with residuals below 



^http: / /heasarc. gsfc.nasa.gov/cgi-bin/Tools/w3nh/w3nh. pi 
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0.6 keV, compared to x^=0.766 for the complex absorber fit. Thus, we prefer the partial 
covering absorber interpretation, also because this is a natural consequence if (as we argue 
below) the X-ray emission region and the absorber are comparable in size: barring a fine- 
tuned geometry, an absorber whose size is comparable to that of the emission region cannot 
completely block our line of sight. 

This intrinsic absorber also offers a partial explanation for the high XSS (>3 keV, un- 
affected by absorption) to RASS (<2 keV, strongly affected) count ratio. However, source 
variability must also play a role. Using the best-fit partial-covering mkcf low model, the 
Suzaku spectrum corresponds to a "predicted" ROSAT PSPC (on-axis) count rate of 0.24 
cts^^. In reality, the vignetting-corrected RASS rate is 0.35 cts~^, so it was likely brighter 
during the RASS epoch by a factor of ~1.5. For the RXTE PCA, predictions (XSS mea- 
surements) are 1.35 cts^^ (2.38) in the 3-8 keV band for a ratio of ~1.7 and 0.51 cts^^ 
(0.96) for ~1.9 in the 8-20 keV band. This could also be due to source variability. Another 
possibility is that there was an additional source within the PCA field of view (but outside 
the Suzaku XIS field of view) during the slew that led to the inclusion of V893 Sco in the 
XSS catalog. 

The 2-10 keV absorbed flux of V893 Sco is 1.7x10"^^ ergscm~^s~^, corresponding to a 
2-10 keV luminosity of 4.9x10^^ ergss^^ at 155 pc assuming isotropic emission. The 0.4-2 
keV absorbed flux is 3.0x10"^^ ergscm~^s~^, while the extrapolated flux of V893 Sco above 
10 keV is < 1 X 10~^^ ergscm~^s~^. This justifies our choice not to analyze the HXD data, 
since a source at this flux level is undetectable given the current systematic uncertainties 
in the background reconstruction. The estimated bolometric X-ray luminosity depends on 
the choice of spectral models (see Tabled]), but may be as high as 1.4x10^^ ergss"^, placing 



V893 Sco among the most X-ray luminous dwarf nova (IBaskill et al. II2005I ). 



3.2. Optical Light Curves 

We had three purposes in mind for obtaining optical photometry of V893 Sco contem- 
poraneously with the Suzaku observations. First, we wished to ascertain the outburst status 
of the object. The instrumental white light magnitudes, as measured with the SAAO 1.9 m 
telescope, were in the range 11.0-11.5, 11.5-12.0, and 11.2-11.8 outside the eclipse, respec- 
tively, on the nights of Aug 25, 26, and 27. Since there is a known offset of approximately 
3.0-3.2 magnitudes with respect to the standard V band for the typical spectrum of a CV, 
V893 Sco had V magnitudes of 14.0-15.0 out of eclipse on these nights. That is, V893 Sco 
was in quiescence during the Suzaku observations. 
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Our second purpose was to c heck the orbit a l ephe meris. We therefore folded the SAAO 
light curves on the ephemeris of iBruch et al. I (120001 ) and found that the eclipse occurred 
~0.17 cycles earlier than predicted. There is a narrow, V-shaped core of the optical eclipse 
lasting ~0.03 cycles, presumably indicating the eclipse of a compact light source such as the 
white dwarf photosphere. We therefore measured the times of the fourth contact, as indicated 
by the change of slope of the light curve, some 0.015 cycle after the mid-eclipse. We also 
measured the times of the first contact, indicated by the slight change in slope during ingress, 
~0.015 cycle before the mid-eclipse, even though there often was ambiguity exactly where 
this happened. The results for the 6 eclipses (2 were observed on each of the three nights) 
are listed in Table[2|, along with the brightness level relative to the faintest point during that 
eclipse. On average, we measured a slight steepening of the gradient at Bruch et al. phase 
0.8093±0.005, ~0.56 mag above mid eclipse; and the end of the steep part of the egress 
at Bruch et al. phase 0.8385±0.0011, ~ 0.47 mag above. W e therefore applied a offset of 
—0.176 to the phases calculated using the lBruch et al. I (120001 ) ephemeris. This is equivalent 
to adopting a linear ephemeris with a period of 0.075961467 d (cf. 0.07596185±0. 00000012 d 
of the ephemeris) while keeping the same epoch for cycle 0, although the time gap of 7 years 
between the Bruch et al. observation and ours is long enough that a quadratic ephemeris 
may be more realistic. 

The final purpose of the optical photometry was a detailed comparison of the X-ray and 
optical light curves. We will present our analysis of the X-ray light curves in §3.3, and report 
on our modeling of both the optical and X-ray light curves in §4.3. 



3.3. X-ray Light Curves 

We present the Suzaku light curve of V893 Sco against both time and relative cycle 
counts (since cycle 34890.0 of Bruch et al. ephemeris, adjusted as described in §3.2) in 
Figure[2l Due to the restrictions of low Earth orbit satellites, the Suzaku coverage has 
many gaps, resulting unfortunately in little strict overlap with the SAAO photometry. The 
phase coverage is highly uneven: it can be seen that phase interval 0.6-0.9, for example, is 
poorly covered. The Suzaku count rate is variable at all orbital phases, but there are dips 
whenever the phase of the optical echpse is covered. The Suzaku data cover the eclipse phase 
completely twice (at relative cycle count 3.0 and 4.0); in addition, they cover a part of the 
ingress once (5.0, as the satellite was approaching the night-side limb of the Earth), and the 
egress once (8.0, as the satellite emerged from the daylight limb of the Earth). 

In the folded representation (Figure[3]), the X-ray light curve shows a significant dip at 
the phase of the optical eclipse. It is also clearly variable at other phases, but no other dips 
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of comparable significance are seen. Given the poor phase coverage of the Suzaku data, it 
is premature to consider any orbital modulations outside the eclipse phase. The apparently 
greater amplitude of variability during the second half of the orbital cycle (0.5-1.0), in 
particular, is probably due to the poorer coverage of these phases. Concentrating on the dip 
at phase 0.0, we investigated its properties in more details to ascertain if it is an eclipse. 

First, we investigated if this dip is energy dependent. In the left half of Figure HI we 
show the folded light curve near the eclipse phase in three energy bands, as histograms. 
Solid lines, repeated in all three panels, are the eclipse profile for the total band. Both the 
individual and the total-band light curves are simply scaled so that the top of the panel 
corresponds to 1.2 times the maximum value of all phase bins for that curve. It is clear 
from these panels that there is little, if any, energy dependence in the depth of this dip. We 
therefore conclude that this is not an absorption event, e.g., by the outer rim of the disk, or 
by an extended atmosphere of the secondary. 

We also investigated how repeatable the dip was from cycle to cycle. A comparison 
between the average profile (shown in Figure[2]) and individual cycles would be dominated by 
the slow (> 1,000 s) variations in the source brightness. For each segment around the eclipse 
phase, we first calculated the average intensities during phases 0.96-0.98 and 0.02-0.04. We 
then normalized the individual light curves using the ratio of these "out-of-eclipse" averages, 
relative to cycle 4. We also calculated the average of the light curves thus normalized. In 
the right half of Figure HI we show both the normalized individual profiles (black histogram 
with error bars) and their average (red line). We confirm that the reduction in count rate 
is seen every time the eclipse phase is observed. At the same time, there are variations in 
the light curves from cycle to cycle. This is most apparent in the egress for cycle 4, for 
which two interpretations are possible. Either the egress was delayed by ~0.008 cycle or the 
total emission was lower by 30-40% during the egress. Given the small number of cycles 
observed, we must keep in mind the possibility that the ingress is similarly affected by a 
similar cycle-to-cycle variability. 

To summarize, we observe a significant decline in the X-ray count rates during the 
phase of the optical eclipse (Figure[3]). This feature is present every cycle, even though 
there are variations in the detailed profile from cycle to cycle, and the fractional decrease is 
independent of the photon energy (FigureH]). We conclude that we have discovered a partial 
eclipse of the X-ray emitting region by the secondary in V893 Sco. 
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Discussion 



4.1. Intrinsic Absorber, Abundances, and Luminosity 



The X-ray spectrum of V893 Sco indicates the presence of a partial-covering intrinsic 
absorber (Tabled]) with Nh ~ 2 x 10^^ cm~^. Such an absorber reduces the observed count 
rate in the ROSAT band, but ha s little effect above 2 keV. Our finding is consistent with 
the study of Ivan Teeseling et al. I (119961 ). who find that intrinsic absorption in the ROSAT 
band is a general phenomenon among high inclination non-magnetic CVs. They conclude 
that this absorber is likely located in the inner disk near the boundary layer. Since the 
innermost region of the accretion disk has a similar physical dimension as the boundary 
layer, partial-covering absorption is a likely result. 

This interpretation is i n conflict with the predictions of the basic version of the disk 
instability model (see, e.g., iLasota 1 1200 ll ) that typically predicts a surface density of the 



inner disk of order 1 gcm~^. Viewing the putative boundary layer through such an inner 
disk should results in Nh ~ 10^^ cm~^. If this is true, the inner disk is optically thick 
to >7 keV photons, and will produce reflection signatures, including a strong 6.4 keV line 
(jPone fc Osborne 1119971 ). The weakness of the observed 6.4 keV line in V893 Sco excludes 



the possibility that such a dense (~1 gcm~^) disk subtends a large (~ 27t steradian) solid 
angle. Our discovery of a ~ 2 x 10^^ cm~^ partial covering absorber in V893 Sco suggests 
that the disk may be far less dense than the basic version of the disk instability theory 
predicts. Note that this same model is already in conflict with ot her key observ ations, such 
as the interoutburst intervals and the quiescent X-ray luminosity (ILasota Il200ll ). 



We note, in passing, that iRamsay et al. I (120081 ) classified three CVs as likely IPs based 
largely on the complex absorber found in their X- ray spectra. Based on our result on 
V893 Sco, as well as the similar finding on OY Car (jPandel et al. II2005I ). we urge caution 
against using X-ray absorption alone to claim an IP classification. 

We infer that the bolometric X-ray luminosity of V893 Sco to be 1.0-1.4x10^^ ergss"^ 
at the time of the Suzaku observations. V893 Sco therefore appears to be among the most 
X-ray luminous dwarf novae. Byckling et al. (in preparation) discuss the X-ray luminosity 
function and temperatures of nearby dwarf novae, including V893 Sco. 



4.2. The Partial X-ray Eclipse 



The detection of a partial eclipse in a non-magnetic CV is unprecedented. All previous 
observations of deeply eclipsing, quiescent dwarf novae showed flat-bottomed X-ray eclipses 



- 12 - 



the corona of the secondary ( 


Wood et al. 1 


1995: 


van Teeseline 


Pratt et al. 


1999: 


Ramsav et al. 


2001; 


Wheatlev k West 


2003) 



1997 



out of eclipses have been used to constrain the total size of the X-ray emission region to 
be not much larger than the white dwarf. On the other hand, the partial X -ray eclipse of 
y893 Sco is also in ma rked contrast to those of U Gem (ISzkody et al. Ill996l ) and WZ Sge 
( jPatterson et al. Ill998l ). These are both lower inclination systems that exhibit a partial 
eclipse of the outer accretion disk and the bright spot in the optical. There is no X-ray 
eclipse in these systems, instead showing energy-dependent absorption dips prior to the 
eclipse. These dips are believed to be caused by photoelectric absorption in an azimuthal 
structures on the outer accretion disk. 

X-ray emission from regions other than the vicinity of the white dwarf and from the 
corona of the secondary have been considered in the literature. One speculated source is 



the br ight spot region. However, this possibility was comprehensively refuted by iPringle 



( 119771 ): in particular, the highest shock temperature possible for the bright spot is of order 
kT~l keV, and therefore this can be discarded for V893 Sco. The other source is the soft 
X-ray component observed in high accretion rate systems, whose physical origin is unclear. 
For example , in the XMM-Newton observations of the eclipsing nova-like system, UX UMa, 



Pratt et al. I (120041 ) found two components of X-ray emission. One is highly absorbed (Nh 
~ 6 X 10^^ cm~^) thermal (a single temperature mekal fit gives kT=5.5 keV) component 
and undergoes a total eclipse, and therefore originates from a compact region around the 
white dwarf. The other component that dominates below 2 keV unabsorbed and does not 
show total or partial eclipse, and is therefore from a highly extended region. However, 
this component has the wrong spectral shape and its origin is far too extended to offer 
an explanation for the partially eclipsed component in V893 Sco. We therefore seek an 
explanation of the partial X-ray eclipse assuming a compact X-ray emitting region centered 
on the white dwarf. 

We believe that a grazing geometry in which the entire photosphere of the white dwarf 
is never completely eclipsed can explain both the optical and the X-ray eclipse profiles. 
To demonstrate that such a geometry is not too contrived, let us consider a plausible set 
of system parameters and estimate the probability of a partial eclipse of the white dwarf. 
Given the 1.82 hr orbital perio d, an assumed m ass ratio of 0.25, an empirical mass-radius 
relati onship of the secondarv (jPatterson I Il984j ) , and a white dwarf mass-radius relation- 



ship ( jPringle fc Webbink I Il975l ). we estimate a secondary radius R2 of 1.28x10^° cm, a 



primary radius Ri of 8.86x10^ cm, and a binary separation a of 4.73x10^° cm (see TableE]). 
With these representative numbers, and approximating the Roche-lobe filling secondary as 
a sphere, the bottom of the white dwarf is just eclipsed at an inclination angle i of 73.88°; 
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the top at i=75.88°, compared to the range of 75.88°-90° for a total echpse. The key factor 
is the ratio of the white dwarf diameter to the radius of the secondary, which is ~l/7 in this 
case, and this should be the approximate ratio of systems with partial and total eclipses of 
the white dwarf photosphere. 

All available data are consistent with our interpretation that the white dwarf in V893 Sco 
is not totally eclipsed. A total eclipse would result in a flat-bottomed optical eclipse. Systems 
with flat-bottomed light curves are often indicated in catalogs as having a double echpse: 
this is the terminology for eclipses that show two episodes of rapid ingress and egress, one 
caused by the total eclipse of the white dwarf and the second by that of the bright spot, 
offset in phase because the bright spot is located away from the line of centers of the binary. 
The white dwarfs in double-eclipsing systems must undergo a total eclipse, because otherwise 
it would be impossible to separate the ingress of the white dwarf from that of the bright 
spot. On the other hand, the optical eclipse in V893 Sco is neither fiat-bottomed or double 



( iBruch et al. 1120001 : this work). Therefore, if the white dwarf is eclipsed at all, it must be a 



partial eclipse. 

We now consider the X-ray eclipse in such a system assuming three different geometries 
of the X-ray emission region: Magnetic pole region(s) as appropriate for IPs, a spherical 
corona that covers the entire white dwarf photosphere, and an equatorial boundary layer. 
We sketch a possible geometry for the boundary layer case in FigureO 

The accretion in IPs, by definition, o ccurs onto an a rea smaller than the entire surface 



of the white dwarf. For the IP, XY Ari, iHellier I (119971 ) used the rapidness of the eclipse 



ingress (<2 s) and their phase jitters (individual orbital phases of the third contact varied 
in the range 0.0543-0.0553) to infer small (<0.002 of the total surface area) spots that move 
on the projected face of the white dwarf as a function of the spin phase. It's extremely 
unlikely that such a small spot be partially eclipsed by the secondary. If, on the other hand, 
two larger spots are located near the rotational poles, a fiat-botto med partial eclipse is a 



likely result, as we do indeed observe in EX Hya (iRosen et al. Ill99ll ). Based on the observed 
X-ray eclipse profile, and the lack of phase jitters, we believe that an IP or IP-like geometry 
is highly unlikely for V893 Sco. 

Next, we consider a spherical corona surrounding the white dwarf. This simple scenario 
has the advantage that, if the white dwarf photosphere is partially eclipsed, so must be the 
corona. Because such a corona is bigger then the white dwarf in every direction, the X-ray 
eclipse must last longer. The depth of the X-ray eclipse should be similar to that of the 
white dwarf photosphere, if the latter can be reliably extracted from the optical eclipse. 
Finally, the X-ray eclipse should be U-shaped, from the limb brightening of an optically thin 
emission plus the fact that the body of the white dwarf blocks the emission from the far side. 
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Although none of these predicted indicators is seen in V893 Sco, the quahty of the current 
data quahty cannot rule out a spherical corona model. 

Finally, we consider the boundary layer picture. We must first consider whether a partial 
eclipse of an equatorial boundary layer is likely. This turns out not to be a severe problem. 
The reason for this is that, due to the inclination of ~ 75°, the projected equator of the 
white dwarf is significantly curved, spanning ~l/8 of the white dwarf diameter. Thus, even 
a smallest possible boundary layer of a negligible height and a negligible extent in latitude 
will suffer a partial eclipse in ~l/8 of the partial eclipsing dwarf nova. The probability can 
only rise with any spatial extent of the boundary layer. 

As a pilot study to demonstrate that a partial eclipse interpretation is valid and that 
future investigations hold the potential to constrain the shape and the size of the X-ray 
emission region, we have created a numerical model of the optical and X-ray eclipse. We 
represent the white dwarf photosphere, a boundary layer, or a spherical corona using a set of 
grid points. We trace the line-of-sight from each point at various orbital phases to determine 
whether it is visible from Earth, blocked by the body of the white dwarf, or eclipsed by 
the secondary. We use the same code to produce the schematic diagram (Figure^]) and the 
model hght curves (Figure[6]). 

For the optical eclipse, we plot the folded optical light curve, normalized to 1.0 at the 
phase of the first contact of the compact light source (at Bruch et al. phase 0.8093; Table[2]). 
In this initial attempt, we only model the eclipse of the white dwarf photosphere, since there 
is little information on the location of the bright spot to allow us to construct a reliable 
model. The eclipse light curve of the photosphere is estimated at 2,000 phase bins. The 
Earth facing side of the photosphere is divided into 90 zones from the center to the limb, 
and 360 points in each zone. The parameters of the light curve model are the Ri/a, R2/a, 
and the inclination angle i. These parameters can be narrowed down using astrophysical 
co nsideration s . For the secondary, we have assumed the empirical mass-radius relationship 



of [Patterson I (119841 ) as before (Table[3]). We then assumed a mass ratio q = M2/M1 = 0.25, 
which determines the white dwarf mass, and hence the radius. We then stepped through 
a range of i to reproduce the phases of the first and fourth contacts (Table[2]). We find an 
inclination of i=74.2° under these assumptions. The model shown reflects the additional 
assumption that the white dwarf was the dominant light source at first and forth contacts, 
and that other sources of light are either too faint (likely to be the case for the secondary and 
possibly also for the accretion disk) or already fully echpsed. Different conclusions would 
result if we assume a different value of q, or if the narrow eclipse is a combination of white 
dwarf and bright spot eclipse. As for limb darkening, we have also s imulated a case with a n 



extreme limb darkening {u=l following, and using the notation of, IWood &: Horne Ill990l ). 
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but this does not make a large enough difference to matter here. The vertical dashed lines 
in Figure[6] are the ingress and egress timings we measured and used to adjust the orbital 
ephemeris (§3.2 and Table[2]). 

For the X-ray light curve, we use the same normalized and averaged version as shown in 
FigureHl For the spherical corona model (shown in green), the parameters are Ri/a, R2/a, i, 
and the height of the corona above the white dwarf h. The corona is divided into 10 layers in 
height, 180 degrees of latitude, and 360 degrees of longitude, and assumed to emit uniformly 
throughout this volume. The system parameter values are the same as those used for the 
optical eclipse, and we assume a height of 0.2 Ri. Changing the height will only affect the 
width of the X-ray eclipse under the assumption of a spherical corona. 

For the boundary layer, we assumed a boundary layer that extends ±5° from the equator 
in latitude (divided into 21 zones), and an X-ray emission region height of 0.2Ri (10 layers), 
which resulted in the red curve. Both the height and the latitudinal extent affects the depth 
and the width of the predicted eclipse, but in different directions in the width-depth plane. 

These preliminary models reproduce the duration and the depth of both the optical 
and X-ray eclipses to better than a factor of 2. To proceed to the next level of actually 
fitting the optical and X-ray light curves, several things must happen. First, the system 
parameters must be constrained from UV, optical, and IR observations to narrow down the 
parameter space. The mass ratio g is a crucial parameter: once we know q, the eclipse light 
curves will in principle allow us to tightly constrain i. It also provides a tight constraint 
on the white dwarf mass, and hence its radius. Knowing q also fixes the initial trajectory 
of matter leaving the Li point, hence constrains the possible location of the bright spot. 
Alternatively, UV sp ectroscopy has the potential to constrain the white dwarf parameters 



( iGodon et al. II2008I ). UV photometry may better isolate the white dwarf eclipse from the 
bright spot eclipse than optical photometry, and allow us to confirm or refute the partial 
eclipse of its photosphere. 

Second, the quality of the X-ray light curve must improve. The statistical quality is 
one issue, but we also need to observe an increased number of X-ray eclipses to average out 
cycle-to-cycle variations. Observations from high Earth orbit will help in allowing a better 
modeling of the longer-term (minutes to hours) variability and to separate that from the 
eclipse profile. 
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5. Conclusions 

We performed Suzaku observations of the eclipsing dwarf nova, V893 Sco. We confirm 
that it is an X-ray luminous dwarf nova, and moreover, report the discovery of a partial 
X-ray eclipse. 

In the past, both optical and X-ray observers have concentrated on systems that exhibit 
a total eclipse of the white dwarf. They provide a strong constraint on the height of the X-ray 
emission region above the white dwarf. However, eclipse light curves are one-dimensional, 
and one can only derive one-dimensional constraints from the eclipse analysis. The X-ray 
emission region, on the other hand, has two dimensions (height and latitudinal extent) 
assuming azimuthal symmetry. It is therefore necessary to observe an ensemble of eclipsing 
dwarf novae, in which the limb of the secondary cuts across the emission region from different 
angles, to be able to constrain the height and the latitudinal extent simultaneously. Inclusion 
of a partially eclipsing system in such a study will be a huge step forward in our quest to 
constrain the X-ray emission region size in two dimensions. Future X-ray observations of 
V893 Sco, together with improved estimates of system parameters from other wavelengths, 
have the potential to test the boundary layer picture of X-ray emission in quiescent dwarf 
novae, and furthermore to constrain the latitudinal extent of such a boundary layer. 

This research has made use of data obtained from the Suzaku satellite, a collaborative 
mission between the space agencies of Japan (JAXA) and the USA (NASA). 
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Fig. 1. — The Suzaku XIS spectra (blue: sum of three units with FI CCDs; red: XISl data, 
containing the BI CCD) of V893 Sco fitted with the coohng flow (mkcf low) model. The 
data are shown in the top panel. The middle panel shows the residuals (as as) when only a 
simple absorber is used. The bottom panel shows the residuals with a simple and a partial 
covering absorber. The solid line on the top panel shows the latter model. The inset in the 
top panel gives an expanded view of the Fe Ka complex, here fitted with the mkcf low model 
plus a 6.4 keV Gaussian. 
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Fig. 2. — The X-ray light curve (black) of V893 Sco in 16 s bins plotted against time (bottom 
axis) and, equivalently, relative cycle counts (top axis). The Aug 26 SAAO photometry is 
shown in blue. Red tick marks shows the two full and two partial coverage of the eclipse 
phase by the Suzaku XIS. 
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Fig. 3. — Folded light curves of V893 Sco from our SAAO (top) and Suzaku (bottom) 
observations. Two complete cycles are shown for clarity. Optical data have been converted 
to linear intensity scale, with a normalization such that 1.0 roughly corresponds to V=14.5, 
or ~125,000 counts in a 6 s integration. 
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Fig. 4. — (Left) The folded X-ray light curves near the eclipse phase in three energy ranges. 
In each panel, the black histograms show the data in the energy range shown, while the red 
line is the total (0.3-10 keV) light curve. All curves are scaled such that the top of the panel 
is 1.2 times the largest bin. (Right) The individual X-ray light curves around the eclipse 
phase. In each panel, the black histograms show an individual cycle (cycles 5 and 8 are 
shown together in the bottom panel). Those of cycles 3, 5 and 8 are scaled by the factors 
shown and explained in the text. Red solid line, repeated in all 3 panels, is the average of 
these normalized light curves. 
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Fig. 5 



— A schematic diagram of the possible geometry of the grazing echpse in V893 Sco. 
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Fig. 6. — Optical (top) and X-ray (bottom) eclipse light curves plotted with representative 
models. The folded optical light curves are normalized at the estimated time of the first 
contact of the white dwarf. The X-ray light curve is normalized so that the average rate at 
phases 0.96-0.98 and 0.02-0.04 is 1.0. For the details of the model parameters, see text. 



Table 1. Spectral Parameters of V893 Sco. 



Model 








CF^ 


kT/kT^,, Abund.^ L/ 








(1022 cm-2) 




(keV) (ergs s~^ 



mekal 1.349 9.34+°i? 0.87+^-^^ l.OTxlO^^ 



mekal 0.951 2.86l°|3 1.96t[J|t 0.51^°;°^ Q.gbt'oil 0.611°;°^ 1.15x10^2 
mkcfiow 1.349 QAQ^oaI 32.4+2-9 l.OStom 1-42x1032 

mkcfiow 0.928 2.68+°-29 1.85l°?^ 0.60+°'°^ 14.5 j:}^ 0.731°°^ 1.28x10^2 



'^The equivalent hydrogen column density of the simple absorber 

^'The column density and the covering fraction of the partial covering absorber 



The global abundance of the plasma, assuming the Solar values of lAnders fc Grevesse 



([19891) 



"^Extrapolated X-ray luminosity in the 0.01-100 keV range, which is a good approximation 
of the bolometric X-ray luminosity, assuming a distance of 155 pc. 
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Table 2. Optical Eclipse Timings of V893 Sco. 



Eclipse Transition 



HJD 



Decline (mag) Bruch et al. Phase 



Aug 


25-1 


In; 


^ress 


2453973 


29557 





63 





81096 


Aug 


25-1 


Eg 


ress 


2453973 


29786 





51 





84113 


Aug 


25-2 


Ini 


^ress 


2453973 


37131 





54 





80805 


Aug 


25-2 


n/g 


i Coo 


9453973 


37367 


n 


41 


n 


8391 3 


Aug 


26-1 


Inj 


gress 


2453974 


28288 





55 





80841 


Aug 


26-1 


Eg 


ress 


2453974 


28531 





64 





84041 


Aug 


26-2 


In; 


gress 


2453974 


35892 





56 





80946 


Aug 


26-2 


Eg 


ress 


2453974 


36122 





45 





83963 


Aug 


27-1 


In; 


^ress 


2453975 


27046 





57 





80936 


Aug 


27-1 


Eg 


ress 


2453975 


27243 





48 





83528 


Aug 


27-2 


In; 


^ress 


2453975 


34644 





40 





80965 


Aug 


27-2 


Eg 


ress 


2453975 


34839 





48 





83525 


Average 


Inj 


2;ress 









56±0.02 





8093± 


Average 


Eg 


ress 









47±0.04 





8385± 



Table 3. Assumed System Parameters 



Mass Ratio"" 


Inclination*^ 


Separation^ Mi^ 


Ri^ 






0.25 


74.2 


4.73x10^° cm 0.58 Mq 


8.86x10* cm 


0.15 Mq 


1.28x10^° 



""Mass ratio is assumed. Inclination is then derived by assuming that the white dwarf 
eclipse has a duration from the first contact to the fourth of 0.0292 cycle. 

'^Separation is calculated using the orbital period and the masses of the two components. 

■^White dwarf mass is derived from the secondary mass and the assumed mass ratio; its 
radius is then derived using the mass-radius relationship. 

"^Secondary mass and ra dii are derived fr om the orbital period assuming the empirical 
mass-radius relationship of [Patterson I (119841 ). 



